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Abstract —In this paper, a new multi-cell MMSE detector is 
proposed for massive MIMO systems. Let K and B denote the 
number of users in each cell and the number of available pilot 
sequences in the network, respectively, with B = /3K, where 
/3 > 1 is called the pilot reuse factor. The novelty of the multi¬ 
cell MMSE detector is that it utilizes all B channel directions 
that can be estimated locally at a base station, so that intra-cell 
interference, parts of the inter-cell interference and the noise 
can all be actively suppressed, while conventional detectors only 
use the K intra-cell channels. Furthermore, in the large-system 
limit, a deterministic equivalent expression of the uplink SINK 
for the proposed multi-cell MMSE is derived. The expression is 
easy to compute and accounts for power control for the pilot 
and payload, imperfect channel estimation and arbitrary pilot 
allocation. Numerical results show that significant sum spectral 
efficiency gains can be obtained by the multi-cell MMSE over 
the conventional single-cell MMSE and the recent multi-cell ZF, 
and the gains become more significant as /? and/or K increases. 
Furthermore, the deterministic equivalent is shown to be very 
accurate even for relatively small system dimensions. 

I. Introduction 

Massive multiple-input-multiple-ouput (MIMO) is an 
emerging technology that scales up multi-user MIMO by 
orders of magnitude compared to current state-of-the-art (T), 
£3- The idea is to employ an array comprising a hundred, 
or more, antennas at each base station (BS) and serve tens 
of users simultaneously in each cell. The system spectral 
efficiency (SE) can be drastically increased without consuming 
extra bandwidth ffl-E). The uplink and downlink transmit 
power can also be reduced by an order of magnitude since 
the phase-coherent processing provides a comparable array 
gain 0. In the limit of an infinite number of antennas, intra¬ 
cell interference and noise can be averaged out by simple 
coherent linear transceivers, and the only performance limita¬ 
tion is pilot contamination and distortion noise from hardware 
impairments 0, 0. These features make massive MIMO 
one of the key techniques for the next generation wireless 
communication systems. 
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In the uplink reception, the most commonly used linear 
detection schemes are matched filtering (MF), zero forcing 
(ZF) and minimum mean square error (MMSE). We consider a 
multi-cell network with B available orthogonal pilot sequences 
and K users in each cell, where B = (3K with /3 > 1 being the 
pilot reuse factor (i.e., only 1//3 of the cells use the same set of 
pilots). In conventional massive MIMO, the BS first listens to 
the uplink pilot transmission from its own users, estimates the 
channels and then constructs user-specific detectors based on 
the channel estimates 0-0. However, more channel state 
information (CSI) can be extracted when /3 > 1. If the BS 
is aware of all pilot sequences, then it can locally estimate 
B channel directions by listening to the pilot signalling from 
all cells instead of only from its own cell. In principle, the 
BS is then able to select its detectors to suppress parts of 
inter-cell interference, since its K users only occupy K out 
of the B channel directions. Based on similar observations, 
the authors of ED propose a multi-cell ZF detector (referred 
to as full-pilot ZF detector in their paper) which exploits and 
orthogonalizes all available directions to mitigate parts of the 
inter-cell interference. The detector in ifTOl achieves higher SE 
than the conventional ZF when the interfering users are near 
to the intended users of a cell. In general cellular networks, 
however, the gain is less obvious, partly due to a loss in array 
gain of B in multi-cell ZF, instead of K with conventional ZF. 
Multi-cell MMSE detectors are proposed in (9) and ifTTl . but 
the former is limited to /3 = 1 and equal power allocation and 
the latter is based on the unrealistic assumption that perfect 
CSI is known at each BS. There is thus need for a more 
practical multi-cell MMSE detector. 

In this paper, we propose a new multi-cell MMSE detector 
which accounts for power control for the pilot and payload, 
imperfect channel estimation and arbitrary pilot allocation. 
By utilizing all the estimated channel directions at a BS, 
the proposed multi-cell MMSE detector can actively suppress 
intra-cell interference, parts of the inter-cell interference and 
noise. Numerical results show that significant sum SE gains 
can be obtained by multi-cell MMSE over the single-cell 
MMSE and the multi-cell ZF from 051 • Furthermore, a 
deterministic equivalent expression of the stochastic uplink 
SINR is derived in the large-system limit. The expression is 
easy to compute and only depends on large-scale fading, power 


control and pilot allocation. The expression is shown to be a 
very accurate approximation even for relatively small system 
dimensions. 

Notations : Boldface lower and upper case symbols repre¬ 
sent vectors and matrices, respectively. The trace, transpose, 
conjugate, Hermitian transpose and matrix inverse operators 
are denoted by tr (-), (-) T , (•)*, (•) H and (•) 1 , respectively. 


K. It is shown later that a larger /? brings a lower level of 
pilot contamination, since a smaller fraction of the cells use 
the same pilot sequences as the target cell. 

Based on the received signal in Eqn. 0, the MMSE 
estimate of the uplink channel h jik is lfl2l 

h jik = Vrndj (zjfc) Yj V* fc , (4) 


II. System Model and Detector Design 

We consider a synchronous massive MIMO cellular network 
with multiple cells. Each cell is assigned an index in the cell 
set C, and the cardinality \C\ is the number of cells. The BS in 
each cell is equipped with an antenna array of M antennas and 
serves K single-antenna users within each coherence block. 
Assume that this time-frequency block consists of T c seconds 
and W c Hz, such that T c is smaller than the coherence time of 
all users and W c is smaller than the coherence bandwidth of all 
users. This leaves room for S = T c x W c transmission symbols 
per block, and the channels of all users remain constant within 
each block. Let h jik denote the uplink channel response from 
user k in cell l to BS j within a block, and assume that it is 
a realization from a zero-mean circularly symmetric complex 
Gaussian distribution: 


hjjfc ~ CAf(Q, dj (zik) Im) • (1) 

The vector z u~ G M 2 is the geographical position of user k 
in cell l and dj (z) accounts for the channel attenuation (e.g., 
path loss and shadowing) related to any user position z G R 2 . 
Since the user position changes relatively slowly, dj (zik) is 
assumed to be known at BS j for all l and all k. In what 
follows, the uplink channel estimation is first discussed to lay 
a foundation for the novel detector design. 


A. Uplink Channel Estimation 

In the channel estimation phase, the collective received 
signal at BS j is denoted as Yj G C MxB where B is the 
length of a pilot sequence (B also equals to the number 
of orthogonal pilot sequences available for the network, as 
mentioned in Section U). Then Yj can be expressed as 

K 

Yj = EZ y/pikhjikv ilk + Nj, (2) 

ze£ fc=l 

where h jik is the uplink channel defined in 0} and pik is the 
power coefficient for the pilot of user k in cell l. The matrix 
Nj contains independent elements which follow a complex 
Gaussian distribution with zero mean and variance a 2 . We 
assume that all pilot sequences originate from a predefined 
orthogonal pilot book, defined as V = {vi, ...,vb}, where 


v 6^ v 6 2 


B, h\ — b 2 , 
0, h\ ^ b 2 , 


(3) 


where is the covariance matrix of the vectorized received 
signal vec (Yj) and is given by 

K 

^ = ( z « m ) v W v iL + v 2i b- (5) 

m- 1 


Then according to the orthogonality principle of MMSE 
estimation, the covariance matrix of the estimation error 

h jik = h jik - hji k is given by 

C jik = E jhjj fc hjifcj 

= dj (z i k ) (1 - pi k dj (z ik) ’I r 7 1 v iife ) I M - (6) 


Notice that 


v"*- 1 = 

Hk 3 


Em=l P(mdj (Zf m )v. 


ilk 


-v H 

9 'I'lk 


r H 


(7) 


where aji is defined to be used later on. Then the estimation 
error covariance matrix in © can also be expressed as 


^jik — dj (z/fc) (1 pikdj (z^) oiji lk B^j Im■ (8) 

As pointed out in JTO], the part Y y (ll> j ) V L of the 
MMSE estimator expression in 0 depends only on which 
pilot sequence 1 ,/k that the particular user uses. Therefore, the 
estimated channels of users who use the same pilot will have 
the same direction, while only the amplitudes are different. To 
show this explicitly, define the M x B matrix: 


H v ,j = 


Wji, 


IVjB 


= Y ,(**) 1 K,...,v|], (9) 


which allows the channel estimate in 0 to be reformulated 
as 

hjifc = pikdj {zik) Hy ; je^ fc , (10) 


where e, denotes the ith column of the identity matrix 1^. 
The fact that users with the same pilot have parallel estimated 
channels is utilized to derive a new deterministic equivalent 
expression of the stochastic SINR in the sequel. 

Notice that the estimated channel h jik is also a zero-mean 
Gaussian vector, and its covariance matrix &jik G (C MxM is 


and define iik G {1,..., B} as the index of the pilot sequence 
used by user k in cell l. In our work, arbitrary pilot reuse 
is supported by denoting the relation between B and K by 
B = /3K, where /? > 1 is the pilot reuse factor and S > B > 


&jik = dj (z i k ) Im — C jik = pikd 2 j (zik) atji lk BI M - (11) 

Define the covariance matrix of hy,ji as &v,ji G C MxM . 
Then according to (ITOl) and <TTTb . &v,ji = ay, HIm. 
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B. Uplink Multi-cell MMSE detector 

After the channel estimation, each received signal at BS j 
during the uplink payload data transmission phase is 

K 

= E E Vn^jlkXik + n j, (12) 

l£C fc=1 

where xik ~ CAf (0, 1) is the transmitted signal from a 
Gaussian codebook and n j ~ CAf(0, ct 2 Im) is additive white 
Gaussian noise (AWGN). Ti k is the transmit power of the 
payload data from user k in cell /, and we use different 
symbols for the pilot and the payload data to allow for different 
power control policies for them. Denote the linear detector 
used by BS j for an arbitrary user k in its cell as g jk, then 
the estimate Xp- of the signal Xjk is 

K 

Xjk = sfkYj = g fk EE VnkiijikXik + gfk^j- (13) 

lec k =i 

Then the following uplink ergodic SE can be achieved 12) 

Rfk = (i - f) E{£ 0)} {log 2 (1 + vfk) } , (14) 

where the SINR rjp. is given in (fl5l ) on top of this page 
and E| h , is the expectation with respect to the channel 

estimates known at BS j. Similarly, E{■|hf J j} in (fl5l> denotes 
the conditional expectation given the estimated channels at BS 
j. The rate in (fl4l) is achieved by treating gj|h jjk as the true 
channel in the decoding, and treating interference and channel 
uncertainty as worst-case uncorrelated Gaussian noise; see ED 
for further details. Thus, R^j, is a lower bound on the uplink 
ergodic capacity. 

The second line of Eqn. (fl~5l> shows that the uplink SINR 
takes the form of a generalized Rayleigh quotient. Therefore, 
a new multi-cell MMSE (M-MMSE) detector that maximizes 
the SINR in (fl5l > for given channel estimates is derived as: 

M-MMSE 

g jk 

= ( E E Tlk ( h Afch)lfc + C jlkj + v 2 Im j h jjk- (16) 

\iec fe=l / 

Furthermore, this detector minimizes the mean square error in 
estimating Xjk as well H1 31 : 

Ej|%fe - Xjk\ 2 |h (i) | . (17) 


By plugging ® and (IT0l> into ( fl6l >. the M-MMSE detector can 
also be expressed as 

g M-MMSE = (,l v , A; n'! ( + (o 2 + <pj) IM)~\jk, 

(18) 

K 

where A j = J2 J2 Ti k pi k d 2 (z ik)e ilk e? k . Apparently, A j is 

/gL fc=1 

a diagonal matrix and the ith diagonal element X t , depends on 
the large scale fading, the pilot power and the payload power 
of users that use the ith pilot sequence in V. The scalar ipj is 


K 

<Pj = EE 7~ikdj (z/fc) (1 pikdj (z/fc) a jilk B), (19) 

IGC k -1 

where aji lk is defined in (|7). 

To elaborate the advantages of our M-MMSE scheme, we 
compare it with related work. First, the conventional single¬ 
cell MMSE (S-MMSE) detector from g)-® is 
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-1 


S-MMSE 

°jk 




I ^ ^ T jm^ 1 jjm^ 1 jjm + Zj + CT Ij^f 


\m= 1 


n jjk! 

( 20 ) 

where inter-cell interference is either ignored by setting Z j = 
0 or only considered statistically as with 


K K 

Z;/ = E < Tj m hjj m hjj m + Tj m hjimhji m > . 

I m=1 l^j m =1 I 

( 21 ) 

Notice that the S-MMSE detector only utilizes the K 
estimated channel directions from within the serving cell, and 
treats directions from other cells as uncorrelated noise. Our M- 
MMSE detector, however, utilizes all the B available estimated 
directions in Hy j so that BS j can actively suppress also parts 
of inter-cell interference when B > K. Therefore, our detector 
can actually maximize the SINR in ( fiTt i. while S-MMSE can 
only do this in single-cell cases. The M-MMSE scheme can 
be seen as a coordinated beamforming scheme, but we stress 
that there is no signaling between the BSs since BS j can 
estimate the channels Hy,j directly from the uplink pilots. 
Moreover, the long-term channel statistics (like the channel 
attenuation and the pilot allocation) of all users to each BS can 
be obtained by using the uplink control channel, and hence 
does not necessarily require BS cooperation. Therefore, the 
M-MMSE scheme is fully scalable. 












Compared with the multi-cell MMSE schemes proposed 
in HD and 0, our M-MMSE detector is more general and 
practical. To begin with, any pilot reuse policy and power 
control are supported in our scheme, which allows for an 
analysis based on a more flexible and practical network 
deployment. As shown in |10], non-universal pilot reuse is a 
reliable way to suppress pilot contamination and achieve high 
spectral efficiency in massive MIMO. In addition, as shown 
later on, it is with the larger pilot reuse factors that the M- 
MMSE achieves large gains over S-MMSE. Furthermore, the 
uplink detector in HD is based on the unrealistic assump¬ 
tion that perfect CSI is known at the BS, while imperfect 
channel estimation is accounted for in our detector. Thus the 
performance gains provided by our detector are achievable in 
practical systems. 

III. Asymptotic Analysis 

In this section, performance analysis is conducted for the 
proposed M-MMSE detector. The SINR in (031 ) is stochastic 
since it depends on the estimated channels in each block. 
Hence, the SE in (11 4b cannot be computed in closed form. 
We compute a deterministic equivalent of the SINR that is 
asymptotically tight. The large-system limit is considered, 
where M and K go to infinity while keeping K/M finite. 
In what follows, the notation M —> oo refers to K, M oo 
such that limsupj U K/M < oo and lim inf mA'/M > 00 
Since B scales with K for a fixed /3, lim sup m B/M < oo 
and liminfjvfS/M > 0 also hold for B. The results should be 
understood in the way that, for each set of system dimension 
parameters M, K and B, we provide a deterministic equiv¬ 
alent expression for the SINR, and the expression is a tight 
approximation as M, I\ and B grow large. The expression 
is “deterministic” because it only depends on the large-scale 
fading, pilot allocation and power control, while the SINR 
in <ED> depends also on the stochastic small-scale fading 
realizations. In what follows, the notation — a ' s - > denotes 

M—>oo 

almost sure convergence of a stochastic sequence. 

Before we continue with our performance analysis, two 
useful results from random matrix theory are first recalled. 
All vectors and matrices should be understood as sequences 
of vectors and matrices of growing dimensions. 

A. Useful theorems 

Theorem 1 (Theorem 1 in IT4lB : Let D £ C MxM be de¬ 
terministic and H £ C MxB be random with independent 
column vectors h b ~ CM (0, Assume that D and the 

matrices R b (6 = 1,..., B), have uniformly bounded spectral 
norms (with respect to M). Then, for any p > 0, 

A„. (d (hh" + - Ltr (DT („)) ^ 0. 

( 22 ) 

lr The limit superior of a sequence x n is defined by limsup n ai n = 
lim (sup : m ^ n}); the limit inferior is defined as liminf n ir n = 

n —>oo 

lim (inf {xm ■ m ^ n}). 

n—>oo 


where T (p) £ C 

T(p) = 


MxM is defined as 

_ f i 

■ + (p) 



( 23 ) 


and the elements of 8 (p) = [<5i (p ),..., 5b ( p)] T are defined 
as 5 b (p) = lim^oo 5^ (p), 6 = 1,..., B, where 


for t = 1 , 2 ,..., 


R _RL_ 

{ b \ M h i+ 5 ? _1) 

with initial values d! 0) = 



(24) 

1/p for all b. 


Theorem 2 (see fl4l ) Let 0 £ C MxM be Hermitian non¬ 
negative definite with uniformly bounded spectral norm (with 
respect to M). Under the same conditions for D and H as in 
Theorem Q] 

j^tr (HH fl + P I M )© (hH h + P I M ) 

--J-tr(DT'(p))-p-A0 (25) 

Ad M—t oo 

where T' (p) £ C MxM is defined as 

T' (p) = T (p) ©T (p) + T(p)d-Y' KbSb ^ , T (p). 

Mfri(i + 8 b ( p)f 

(26) 

T (p) and 5 (p) are given by Theorem Q] and 5' ( p ) = 
[<% (p), ..., S' B ( p)] T is calculated as 

S' (p) = (Is-J(p)) _1 v(p) (27) 


where J (p) and v (p) are defined as 

l 


[J (P)] M = 


M tr(R b T(p)R;T(p)) 
M(l+<5 i (p)) 2 


, 1 < 6, l < B 


(28) 


[v (p)] 6 = -jtr (R b T (p) ©T (p)) ,l<b<B. (29) 

B. Deterministic Equivalents of the SINR in ( H5D 

In what follows, we derive the deterministic equivalent fj// 
of rf/l for the M-MMSE detector such that 




rifk 


———— > 0. 
M—t oo 


(30) 


Theorem 3 For the uplink MMSE detector in (fl8t . we have 
> 0, where is given in (31) on top of the 


vfk ~ r lfk 
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where 

1) T j = Tj (a) and S (a) = [<5i,..., 5b} T are given by 
Theorem Q] for a = a ^ Pi and R b = 7 j b $ V j b . 

2) T ' jk = T jk (a) and S'(a) = [5[,S' B f are given by 
Theorem |5] for a = f l f! . © = &v,ji jk , and R/, = 
7 jb®v,jb- 

3) T ” k = T'- k (a) and S' (a) = [5[, are given 

by Theorem [2] for a = 17 ^ Pj , 0 = I m, and R b = 

ljb®V,jb- 

Proof: The main idea behind the proof is to derive the 
deterministic equivalent of each term of the first line of (IT5l) . 
Then the deterministic equivalent of (IT5ll is obtained as (31). 
The full proof can be found in the Appendix B of ED. ■ 
With the above, the ergodic SE in (ITdl i. after dropping 
the prelog factor, converges to Rj k = log 2 (l +vjl)- Since fj^ k 
does not depend on the instantaneous small-scale channel fad¬ 
ing, (1 — §)Rf k is a large-scale approximation of the ergodic 
SE. As shown later, this approximation is even very accurate 
at small system dimensions. Furthermore, the approximation is 
easy to compute and allows for simple performance analysis. 

IV. Simulation Results 

In this section, we illustrate the accuracy and usefulness of 
the analytical contributions for a symmetric hexagonal network 
topology. We apply the classic 19-cell-wrap-around structure 
to avoid edge effects and guarantee the same performance for 
all cells. Each hexagonal cell has a radius of r = 500 meters, 
and is surrounded by 6 interfering cells in the first tier, and 
12 in the second tier. To achieve a symmetric pilot allocation 
network, the pilot reuse factor can be /3 £ {1, 3,4, 7,...}. 

User locations are generated independently and randomly in 
the cells by following uniform distributions, but the distance 
between each user and its serving BS is at least 0.14r. For each 
user location z £ R, a classic pathloss model is considered, 
where the variance of channel attenuation is d , (z) = -m— ( , . 
Here b, £ M 2 is the location of the BS in cell j , k is the 
pathloss exponent, and ||-|| denotes the Euclidean norm. C > 0 
is independent shadow fading with 10 log 10 (C) ~ Af( 0 , cr 2 ^). 
We assume k = 3.7, a 2 / = 5 and coherence block length 


S = 300@ Orthogonal pilots are chosen from a B x B discrete 
Fourier transform unitary matrix. 

Statistical channel inversion power control is applied to both 
pilot and payload data transmission GO), i.e., pi k = rik = 
dl ( , zlk / ) ■ Thus the average effective channel gain between users 
and their serving BSs is constant: E {pi k ||hnfc|| 2 } = Mp, and 
the statistical per antenna received SNR of each user at its 
serving BS is pi a 2 . This is a simple but effective policy 
to avoid near-far blockage and, to some extent, guarantee a 
uniform user experience. In the simulation, p/a 2 is set to 0 
dB to allow for decent channel estimation accuracy. 

In the simulation, 10000 independent Monte-Carlo channel 
realizations for small scale fading are generated to numerically 
calculate the achievable SE in (IT-fl ). The numerical result and 
its large-scale approximation from Theorem 0 are shown in 
Fig. □ for K = 10 and different M. Fig. |T| shows that the 
achievable sum SE increases monotonically as /3 grows, at 
least for /? < 7. This is a result of the following two properties. 
Firstly, a larger /3 results in a lower level of pilot contam¬ 
ination, contributes to a higher channel estimation accuracy, 
and thereby increases the system SE. Secondly, a larger /3 
also indicates more available estimated channel directions in 
the M-MMSE detector, thus a higher inter-cell interference 
suppression can be achieved. Fig.[T|also shows that the Monte- 
Carlo simulations and the large-scale approximations match 
well, even for relatively small M and K. 



Fig. 1. Achievable sum SE as a function of the number of antennas M, for 
(3 e {1,3,4, 7} and K = 10. 

To show explicitly the advantages of our M-MMSE detec¬ 
tion scheme, simulation results for the matched filter (MF) 
from G), the multi-cell ZF (M-ZF) detector from iflO) . and 
the S-MMSE detector from (l20l > are provided for comparison. 
Notice that M — (3K > 0 is needed for the M-ZF scheme, thus 

2 This coherence block can, for example, have the dimensions of T c = 3 ms 
and W c = 100 kHz. 



















the minimum value of M for the M-ZF is /3K+1. Since Fig. |T| 
shows that /3 = 4 and /3 = 7 give the highest performance, 
we provide simulation results for /3 = 4 and (3 = 7 in Fig. [2] 
and Fig. [3] respectively. The MF scheme always achieves the 
lowest performance since it does not actively suppress any 
inter-user interference. Compared with the S-MMSE, our M- 
MMSE achieves a notable SE gain and the advantage becomes 
more significant as 8 and/or K increases. For 8 = A and 
M = 200, the SE of M-MMSE are 28% and 56% higher 
than those of S-MMSE for K = 10 and I\ = 30, respectively. 
When 8 = 7, the gains increase to 40% and 84%, respectively. 
Notice that when 8 = 7, K = 30 brings lower achievable rates 
compared with K = 10, due to the loss from a large pilot 
overhead. The advantage of the M-MMSE over the M-ZF is 
relatively small for K = 10, but it becomes significant as 
8 and K grow. Moreover, the M-ZF can sometimes achieve 
very low SE for small M, while our M-MMSE can always 
achieve good performance, with the same complexity as for 
the M-ZF. From the analysis above, it can be concluded that 
the proposed M-MMSE brings a very promising gain over 
single-cell processing and the M-ZF, and the gain becomes 
increasingly significant as 8 and/or K grow. 



Fig. 2. Achievable sum SE of M-MMSE (squares), M-ZF (triangles), S- 
MMSE (diamonds) and M-MF (circles) with /3 = 4, K = 10 and K = 30. 

V. Conclusions 

In this paper, a multi-cell MMSE detector is proposed and a 
tight deterministic equivalent SINR expression is derived in the 
large-system limit. Compared with the conventional single-cell 
MMSE scheme, that only utilizes the estimated directions from 
within the serving cell, the proposed multi-cell MMSE scheme 
utilizes all channel directions that can be estimated locally 
at the BS to suppress the inter-cell interference. Numerical 
results show that the proposed multi-cell MMSE brings very 
promising sum SE gains over the single-cell MMSE and the 
multi-cell ZF. Since imperfect estimated CSI is accounted for 
in our scheme, the gains obtained are likely to be achievable 
in practical systems. The M-MMSE scheme is the new state- 
of-the-art method for massive MIMO detection and is hard 
to beat since it maximizes the SINR under very general 
conditions. Furthermore, the deterministic equivalent is shown 
to be accurate even for relatively small system dimensions. 



Fig. 3. Achievable sum SE of M-MMSE (squares), M-ZF (triangles), S- 
MMSE (diamonds) and M-MF (circles) with /3 = 7, K — 10 and K = 30. 
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